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Abstract. The Space-Time Integrated Least-Squares (STILS) method is considered to analyze 
a space-time domain decomposition algorithm for scalar conservation laws. Continuous and discrete 
convergence estimates are given. Next using a time-marching finite element formulation, the STILS 
solution and its domain decomposition form are numerically compared. 
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1. Introduction. The need of using domain decomposition methods (DDM) 
when solving partial differential equations is nowadays more and more obvious. Es- 
sentially, it consists in splitting the domain into several parts and in solving a family 
of problems within the subdomains: this family is equivalent to the original problem 
if and only if suitable conditions holds along the interface separating the subdomains. 
This technique is motivated essentially by numerical approximation (for instance, 
computational complexity can be reduced and geometrical features can be exploited 
thoroughly). However, the study of convergence of any discrete scheme usually needs 
a careful investigation of the properties of the continuous solution. The classical way 
of applying domain decomposition methods to evolution problems is to discretize the 
time dimension first uniformly over the whole domain by an implicit scheme, and 
to apply domain decomposition at each time step separately to solve the sequence 
of steady problems obtained from the implicit time discretization. Theoretical and 
numerical analysis of this approach with non-overlapping subdomains for hyperbolic 
problems as transport equation, scalar conservation laws can be found in [551 HZl ES] ■ 
Another way to consider the domain decomposition approach is given by the over- 
lapping Schwarz method which is applied notably to the heat equation see |26) . the 
convection-diffusion equation see [10] and [H] . The possibility to use different grids in 
space for each subdomain is analyzed in [3j. But due to the uniform time discretiza- 
tion, one cannot have an optimal space-time discretization. Thus this approach has a 
significant disadvantage to use a uniform time discretization over the entire domain 
and thus loses one of the main features of domain decomposition methods, namely 
to treat the problem on each subdomain numerically differently, with an appropri- 
ate discretization both in time and space adapted to the subdomain problems. For 
these reasons space-time domain decomposition methods are mandatory, like Schwarz 
Waveform Relaxation (SWR) method, which is widely used. The scheme also uses 
an overlapping domain decomposition in space, like the classical Schwarz algorithm 
for steady problems see [21] , but then the algorithm solves evolution problems on the 
subdomains and uses an iteration to converge to the solution of the original problem. 
Convergence properties of the SWR method applied to convection dominated viscous 
conservation laws with nonlinear flux are given in [I5j . 

This work is a contribution to the non-overlapping case in dealing with a new space- 
time DDM for scalar conservation laws without viscous terms and using the STILS 
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method. The STILS approach for PDEs allows to stabilize the numerical finite ele- 
ment solutions. 

More generally, the method of least-squares is a standard approach to the approxi- 
mate solutions of PDEs see for example in [SUl H] . This method allows to convexify 
noncoercive problems. The basic idea is the same in finite dimensional case: a non- 
degenerate system is transformed in a positive definite system. In |22j . it has allowed 
to perform the Stream Upwind Pctrov Galerkin (SUPG) method [5] in the Galerkin 
Least-Squares (GLS) method. In the same spirit, the Streamline Diffusion and Dis- 
continuous Galerkin [531 HI] have mutated to give Characteristic Streamline Diffusion 
(CSD) [5T]. The space-time framework is presented for example in [T2 HI]. The 
common principle of these methods is to stabilize the problem in adding a quadratic 
residual to the advective terms. Next the STILS method (without to add a quadratic 
term) is introduced in [TJ [31 13] for the transport equation with more convincing re- 
sults. An existence-uniqueness analysis of the STILS solution and his comparison with 
renormalized solutions have been detailed in [S] for linear conservation laws when the 
velocity has a low regularity. Finally in the reference [6] the time-marching version of 
the STILS method has given well numerical results compared to the solutions obtained 
in characteristic method or in the usual finite element methods (Galerkin, Streamline 
Diffusion, Shock Capturing). 

In section 3, the original problem of scalar conservation laws and its two-domain form 
are formulated in a space-time form. These STILS formulations with an equivalence 
result and an iteration-by-subdomain scheme are presented in section 4. In section 5, 
the established results are generalized to the multi-domains case in time. In section 
6, discrete convergence estimates are given with finite element methods in time and 
space. Finally in section 7, the STILS and STILS-DDM solutions are numerically 
compared for the Hansbo example providing in [21j . 

2. Problem statement. Let fl C R*^ be a domain with a Lipschitz boundary 
dfl satisfying the cone property, and let T > 0. Consider an advection velocity 
u : nx]0,T[ — > R'', and / € i^(17x]0,r[) a given source term. In the sequel we 
denote by (.|.) the euclidian scalar product. 

The problem consists in finding a function c : f2x]0,T[— > R satisfying the following 
partial differential equation 

(2.1) g+div(c^)=/ 

and the initial and boundary conditions 

(2.2) c{x, 0) = Co{x) for x in fl 

(2.3) c{x,t) = coo{x,t) for x in dQ- x (0,T) 

where dil^ = {x e dfl : {u{x , t)\n{x)) < 0}, n{x) is the outer normal to dH. at point 
x. For the sake of the presentation, it is assumed that dil^ is not dependent on t. 
The function c represents for example a chemical concentration transported by a liquid 
flow in the domain fl. The velocity u may be a solution of Navier-Stokes equations. 
It may too derive from a potential. For example in a porous medium context, u will 
be a Darcy velocity. 

Let us now partition Q into two non-overlapping subdomains Qi and with Lipschitz 
continuous boundaries dQi and dil2. The common interface dfli ndn2 is denoted by 
F; the normal unit vector on F pointing into is denoted by n. 

(2.4) n^niun2, F = 9r2in9rj2, on- = dfi- ndn„ i = i,2. 
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Then the original problem is rewritten in a two-domains form. More precisely for 
i = 1, 2 we consider the problem 



(2.5) 
with 



dt 



+ div(c'u) = / 



(2.6) c*(a;, 0) = co(a;) for x in fi^ 

(2.7) c\x,t) = coo{x,t) for x on dn^ x (0,r) 

and the following interface condition 

(2.8) = for X onr X (0,T). 

The equivalence between the two above problems is proved in the steady case (see 

mm)- 

The originality of the STILS method is to solve the equation (2.1 ) in both directions 
(space and time) together and in sense. In the next sections the principles of 
the STILS method [D H E E] are recalled for the equations ([2l])-([^ , and its 
two-domains form (|2.5[)-(2.8) is also equivalently formulated. 



3. Space-time framework. In a space-time description, time is just seen as the 
(d-l-l)-th space dimension. Introduce the following space-time objects: the domain 
Q = f]x]0, T[, with the boundary {dnx]0, T[) U (17 x {T}) U (f7 x {0}), and the outer- 
normal n is given by 

on dnx]0,T[ 
on n X {0} 
on nx {T} 

and the space-time velocity is m = (ui, ...,itd, !)*• The space-time inflow boundary is 

dQ~ = {{x,t) e dQ, {u\n) < 0} = (9f7"x]0,r[) U {Q x {0}) . 
In this context, initial condition and inflow condition become 




Cbix,t) 



cq{x) if {x,t) e fl X {0} 
CoQ{x,t) if {x,t) e dfl- X {0,T). 



For a sufficiently regular function ip defined on Q, 

d 



dip 

dx2 ' 



1^,^) IS 

OXd ' Ot ' 



the space-time gradient, and V( • (ipu) 



dip 
dt ' 



d 

■i— 1 



( dip 

{ui(f ). With the above notations. 



the problem (2.1)-(2.3| transforms in space-time form 
(3.1) Vt • (cu) = f 



(3.2) 



Cf, for {x, t) on dQ 



And the space-time counterpart of the two-domains problem (2.5)-(2.8) is given by 
(3.3) Vt • (c'u) = / 



(3.4) 
(3.5) 
(3.6) 
(3.7) 



Cb,l 

/ 

Cb.2 



for {x, t) on dQi 
for {x, t) on r 

for (x, t) on 
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where (see fig. 3.1 ) 

= X (o,r), 



r ^dQiD dQ2 = u r+ u r° 

= {{x,t) e r : {u{x,t)\n{x,t)) < 0} = fj^ = f+ 
f + = {{x,t) e r : {u{x,t)\n{x,t)) > 0} = = f + 
f" = {{x,t) ef : {u{x,t)\n{x,t)) ^0} 




Fig. 3.1. Decomposition of Q and explanation of the symhology 

In the sequel the velocity has the following regularity 

(3.8) ue L°"{Q)\ V-we L°°(Q). 

The following theorem is proved in |13j . 

Theorem 3.1. Under the assumption (3.8) the normal trace of u, (u\n) is in 
L^idQ). 

4. STILS method. 

4.1. Global formulation. For ip E 'D{Q), consider the norm 



m\H{u.Q) = I Wvfmo) + l|Vt • {uLp)\\\2,Q) + / \{u\n)\Lp^dd 

JdQ- 

see [5]. Define the anisotropic Sobolev space H{u,Q) as the closure of I?(Q) for this 
norm 



H{u,Q)^V{Q) 



H(u,Q) 



Under the assumption (3.8), there exists a linear continuous trace operator [5] 



7- : H{u,Q) ^ L^{dQ- Mu\n)\da) 
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Finally define the spaces 

Hoiu, Q, dQ-) = e H{u, Q), = 0} 

V-={^e H{u,Q), 7-(^) e L^idQ-, \{u\n)\dS)} 
G-=j-iV-). 

The fohowmg theorem is proved in [4, 5 . 

Theorem 4.1. (Curved-Poincare inequality) 
If u E L°°{Q)'^ and V • u G L°°{Q), the semi-norm on H{u, Q) defined by 



(Vt • iuip)f dxdt - 



|(w|n)|(y5^d(T 



is a norm equivalent to the norm on H{u,Q). 



A solution of equation ( 2.1 ) in L corresponds to a minimizer in {ip € H{u, Q); 7 (93)- 



Cb = 0} of the following convex, H(u,Q)-coercive functional 



Jic) - \ 



(4.1) 

Its Gateaux derivate is 
(4.2) [DJ{c)]{ip) = 



(Vt • [uc) - fY dxdt 



dQ- 



(?{u\n)(id 



(Vt • {uc) ~f)Vf (uip) dxdt - 



Q 



dQ- 



C(p{u\n)da. 



Therefore a sufficient condition to get the least-squares solution of (2.1|-(2.3| is the 
following weak formulation: if Cf, G , find c G H(u, Q), 7^(c) = such that 

(4.3) Jq Vt • (uc)Vt • {uif) dxdt = Jq^^i- (uip) dxdt Vip G Ho{u, Q, dQ-). 

The maximum principle and weak-solution results proved in 5J are given by 
Theorem 4.2. (Maximum principle) Assume thatfl is bounded, and the function 
/ = m (^.Jp. Let Cb e G- n L°°{dQ~), ifV-u = the solution c of satisfies: 



inf Cb < c < sup Cb- 



Theorem 4.3. (Weak-sohtton) For f G L°^{Q) and u G L^{S),T;BV{n)'^) with 
V • u = 0, let c £ Hq{u, Q, dQ-) be the solution of 

/ Vf • {uc)Vt ■ (uip) dxdt = / /Vt • {u(p) dxdt V(p G Ho{u, Q, dQ-) 

Jq Jq 

then it is a weak solution, i.e: 

/ Vt • (uc)(/3 dxdt = / ftp dxdt V(y3 G i^(g). 
Jq Jq 

For Cb G G- , let Cb G H{u,Q) be such that 7~(Cfc) — Cb. Then p = c — Cb G 
HQ{u,Q,dQ-) is the unique solution (Curved-Poincare inequality (theorem 4.11 and 
Lax-Milgram lemma) of 

(4.4) / Vt • (Sp)Vt • (uf) dxdt = / /Vt • (uv?) dxdt - / Vt • {uCb)Wt ■ {tup) dxdt 
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V(p e Ho{u,Q,dQ^). Set now 

Wiu,Q) = e L^Q), Vt • (uifi) e L^iQ), ipig^_ G L^dQ- , \{u\n)\da)}. 

If u is regular enough, e.g. u £ L^(0, T; i/p (fi)'^), it can be seen that 
HiU, Q) n L°°(Q) = W^(f/, Q) n L°°(Q) 

see |S]. 

In the sequel we assume Cf, € L °°{dQ^) n such that Cf, ~ 7^(Cb) with Cf, G 
Thus, if / = in the least-squares solution c e L°^{Q). 

4.2. Two-domains formulation. In this subsection the two-domains problem 



(jS^-lja/Tj) is rewritten in STILS form. 

4.2.1. Anisotropic Sobolev spaces. Consider the norm 



\Wt\\H(u,Qi) = W'PiWhiQ,) + • {uip^)\\\2^Q^^ + 



\{u\n)\(p^^dcr 



and define H{u, Qt) as the closure of 'D{Q^) for this norm 

Moreover if E is a (relatively) open set of dQ^ U or U , set 



|(u|H)|dCT) = {(^ : S ^ R; / \{u\n)\ip^da < oo} 



this is a Hilbert space with the norm 
induced by the scalar product 



|(u|n)|(/3^diT 



1/2 



For i = 1,2, consider the Hilbert space 

Hidiv, Qi) = {(^, e L'(Q.)''+' : Vt • G i'(Q,)} 
with the graph norm 

1/2 



H(div,Q,) 



ll¥'i|lL2(Q^)d + l + ||Vt • ¥'l||L2(Qi) 



The same definition is given to i?(div, Q). 

The trace result for the H{u, Q^)— functions is given by 
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Proposition 4.4. Under the assumption {S.Sjj there exists a linear continuous 
trace operator 

-ir-H{u,Q,) L^{dQ^,\{u\n)\da) 

which can he localized as 

j,^:H{u,Q,) L^dQfuff,\{u\n)\da) 
'Pi ' — > . ~. 

Proof. The detailed proof can be found in first, consider the well-know trace 
operator 7„ defined from iJ(div, Q^) with value in H~^{dQi) (see [El [7]) 

H> (n|wi)|gQ^, for all Vi e i/(div, Qi) 
with associated Green formula 

/ {Vt■{v^)4>^ + {v^\^ti^^))dl^dl^({v^\n),'^l}^) -i ^ , Vl^, G ^ (Q.) ■ 

Jq. \ I H 2 (dQiy,H2{dQi) 

Plugging Vi = upi in the previous formula, we have 

(Vt • + (^^iVtV'OPi) dxdt = Pj( ) -i i , \/i>i e H'^ {Qi), 

Q. \ / H 2 {dQi);H2 (dQi) 

Let us now consider the bilinear form Li : 'D{Qi)x'D{Qi) C H{u, Qi)xH{u, Qi) — > R 
defined for all ipijipi € 'D{Qi) by 



L,{ip,,'ipi) = I (Vt ■ iuipi)^i + {u\Vt'>JJi)(pi)(ixdt + / _ \{u\n)\(p,->jjtdcr 

JdQT urr 



Accounting to theorem |3.1| we have 

\Li{Vi,^i)\ < ||Vf • (u</?,,)||i2(Q.)||7/>,||i2(Q.) + ||Vt • {uij},) - Vf • {u)i:^\\L2[Q^)\\ipi\\L■^{Q,) 

+ 1 1 1 1 - ur- , h n) 1 1 "^^ 1 1 oq- u?- , (Z\n) ■ 
And the following estimate holds true 

\Li{ipi,i},)\ < (1 + ||V •M|lioo(Q^))||(^,||^(„^Q^)||V'»||/i-(„,Q^). 

Since it 's straigforward to check that Li(ipi,ipi) — , ~, — , if we extend by 

SQ^Ur^ ,(ii|n) 

continuity the bilinear form Li to H(u, Qi) x H{u, Qi) □ 

The Curved-Poincare inequality for the Qi)— functions is given by 
Theorem 4.5. Under the assumption (3.^ the semi-norm on H{u,Qi) defined 

by 



\^i\H{u,Qi) ^ { I [Vt ■ iulp^))^dxdt+ _ \{u\n)\ip'^da 
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is a norm, equivalent to the norm given on H(u,Qi). 

The proof is detailed in [H |S] and is achieved by the inequality ||'y3i||L2(Qi) ^ 
C{T)\(pi\H{u,Qi)^ where C(T) = max{AT^ , 2T) in the divergence free advection ve- 
locity case. 

Consider now the least-squares formulation of the two-domain space-time problem 
([33|-([37l: find {c\c^) G H{u,Qi) x H{u,Q2) such that 



(4.5) 

(4.6) 
(4.7) 



(4.8) 
(4.9) 



/ Vt • (uc^)Vt • {u(p'^) dxdt = / /Vt • {u(p^) dxdt 
JQi Jqi 

= Cb 1 on dQi 

/ Vt • (uc^)Vt • (uip^) dxdt = / /Vt • (u</?^) dxdt 

Jq-2 JQ2 

on 



C = Cb_2 
'r-,(M|n) - 'r+ 

for ah iip\ip^) e Hoiu,Qi,dQ]; uf];) x Ho{u,Q2,dQ^ uf^), and 
for ah^ e A = {/i e L^if- uf+ ,\{u\n)\da), 3 e H{u,Q), 



Ai}nL°°(r- u^^ 



Let us show that (4.5 1-(4.9 1 can be seen as a two-domains form of the STILS problem 

When the flow field has a constant direction at subdomai n iiit e rfac e. for example 
{u\h) is always positive (F^ = 0), the two-domains problem (4.5)-(4.9l has a unique 
solution. Indeed, the equation (4.5)-(4.6) c an b e independently sol ved a nd a dmits 
a unique solution G H{u,Qi) (theorem 4.5). Next one solves (4.7l-(4.8) with 

'r+ 'r+ _ _ 

Let us investigate the situation F" 0, r+ 0. 

It's well known that Green's formula holds for all (pi £ H{div, Qi) and all tpi G H^{Qi) 
(see Lions and Magenes [25]): 



V'i(div Lpi) dxdt 



(V^ilyp,) dxdt 



{lpi\n)'4)^da 



dQi 



Let be a vector valued function of _ff (div, Qi), i — 1,2, and C be the function in Q, 
whose restriction to Qi coincides with Qi- Then C, belongs to 77(div, Q) if and only if 

(Ci|H) = (C2|n) {Hi{f)y. 

Since the STILS solution c e H{u, Q), it's seen that uc e H{div, Q) and ci € i^(F). 

r 

Moreover, we have 

{uc\n) e L^(T- U r+), (Sc^ln) = {uc^\n) (H^oi^- U ?+))'. 

Finally, as Hqq{T^ U F+) C L^{T^ U F+) with dense and continuous injection see 
[ISlIIi], we obtain 



a.e on F~ U F^ 
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and (4.9) is satisfied. 



If c € H{u,Q) = V{Q) 



H{u,Q) 



is the solution of (4.3), then one gets that ci € 

c is a weak solution. 



^—H(u,Q,) 

H{u, Qi) — T^iQi) ■ Under the assumptions of theorem 

and Vt • (5c) — / = in L'^{Q). Let ipi G Ho{u, Qi, dQ~ U F"), one has 



4.3 



I / (Vi • {uc\q) - f)Vt ■ iW) dxdtl < ||V, • {uc) - /||L^(Q)||Vi • {u^')\\l-(q,) = 0. 



Finally, (c|g ,C|q ) is a solution of (4.5)-(4.9). Thus we have the following result. 

Proposition 4.6. Assuming u G 1/^(0, T; BV{^Y) with \/t -u ^ and f e 
L°°{Q), the STILS solution of (4--S) is a solution oj the two-domains STILS problem 
{4-5)-{4-9) i.e: // c G H{u,Q) is the solution of (4--3) then (cIq^iCIq^) G H{u,Qi) x 
ji(u, Q2) satisfies U^-IJ^. 



To get the reciprocal it suffices to show that the interface equation associated to 
(4.9) admits a unique solution. 

In the sequel {^., denotes the duality pairing between W and its topological dual 
W. 



w 



4.2.2. Interface equation. Let 77 G A, and let G Ho{u,Qi,dQ^ ) be the 
STILS-homogeneous extension of 77 on Qi, i = 1,2 defined by 

(4.10) / V( • (5c'''')Vt • (5(^0 dxdt - yip,eHoiu,Q„dQ-uf-) 

JQi 

(4.11) 7r(c'''')|~ ^V- 

^ 

Let also c''* G Ho{u,Qi,T~) be the Dirichlet homogeneous extension satisfying 

(4.12) [ Vt • (uc'^*)Vt • {uip') dxdt = [ f\/f {uLp') dxdt 

(4.13) 7r(c^'*)| =c,. 

'aQ. 

for all (ys' G Hi^{u,Qi,dQ^ U F^). Thanks to Lax-Milgram Lemma, the Curved- 
Poincare inequality (theorem 4.51 shows that these solutions are unique. 
Let (c\c^) a solution of the problem (4.5|-(4.9|. If A = c^|~_^~+ = c^|~_^~+ G A 
then d = c*'^ + d'* . Conversely, {c^'^ + ?'*,^^ + c^^*) G H{u,Qi) x H{u,Q2) is a 
solution of the two-domains formulation if and only if it satisfies (4.9) : 

„1,A 



lAr- rr\ + (A — 



r+,(ti|n) 



= 0, V^i G A. 



Thus to obtain a solution (c , c ) of (4.5)-(4.9) such that G A, i = 1,2 it is 

sufficient to find A G A satisfying the yteklov-Poincare equation 

(4.14) (5A,m)^ = (x,m)^ Va^GA 

with = 5*1 + 5*2 and X Xi + X2 such that 

(5,A,m)^ (A,m)~_^(~|~^ - (c*-\Ai)p+ (~|~), 
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and 



Consider G 'D{Q^) and Q a regular enough function, we have 



Vf • (uCiCi)dxdt = / (u|n)^iCidcr = / _ (M|n)^iCidcr + / _ {u\n)£^iC,ida 
then 

JdQ+ur+ JdQiurr Jq, 



and 



+ / (Vt • (26)0 + / (Vt • (uC.)e< - (Vt • u)e.C.)dxdt 
jQi Jq, 

Thus by density, if Vt • u = the following result is obtained. 

Proposition 4.7. Let G H{u,Qi) such that • (uv'*) = in L'^{Qi), then 



JdQfurf 



If u G L^{0,T,BV{ny), Vt • M = 0, the solution c'''' of the subproblem pOO )- 
(4.11 ) is also a weak solution (theorem 4.3). Thus Vt • {uc^'^) = in L'^{Qi). And the 



above proposition shows that S = Si + S2 ■ A — >■ A' such that Sitj = 77|~_ — c''''|~+ 

is continuous (i.e SiTj G A'). 
Let ?7 G A we have 



I")" 



Moreover 
(4.15) ^11 



T+,(m|ti) 



r+:(«l«) 2 



\\v-c''% 



r+,(M|n) 2 "r+,(«|n) 



Finally applying relation (4.15), and proposition 
result is obtained. 



4.7 



to (^SiT],!]^ , the following 



Proposition 4.8. For 1,2 we have 

^) " '"r-,(n|n) 

Thus for all 77 G A one gets 



\ /a -2^11 "aQ+,(«|n) " '"r-,(n|n) "r+,(n|n) " ' " r+ , (u|n) ^ ' ' 



(^.,.)^>^(I|C^'1^,.,,„) + I 



,2,»?||2 



aQj,(ti|ri) 



r^,(«|«) 



77 — c 



2.»)||2 
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Consider now 77 G A such that Si] ~ 0. Define on Q such that c''|q — d'^ . Then 

to the subproblem 



4.2 



7~(c'' ) = and Vt • (uc'') — in L'^{Q). Applying the theorem 
(4.10)-(4.11) with u regular enough then c** € Ho{u,Q,dQ~) and is the least-squares 
solution of (|3.1[ )-(3.2 ) with / = and Cf, = 0. Using the Curved-Poincare inequality 
(theorem [Zlpwe have c'' in L'^iQ), 77 = G L^j-f- y |(2|n)|d5), so S is 
positive define. Thus the following result is proved. 

Proposition 4.9. Ifu e 2.2(0, T; Vf -u 

ator S realize an isomorphism between A and S{A). 



the Stelov-Poincare oper- 



For / = 0, we have ci 



lr-ur+ 



e A and x G S'(A). Thus the interface equation 



admits a unique solution. 
Finally the following result is obtained. 

Theorem 4.10. Assmning u G L'^{0,T;H^{n)'^), VfU 



Q and J in {2.1) 



the interface equation (4. 14-) associated to admits a unique solution. Thus if 

c€H{u,Q) is the solution of then {c\ , ci ) G (H{u,Qi)x H{u,Q2))^L°°{Q) 

'Wl ' 

satis fi es 14-. ^ -(4.9). Reciprocally if {c^ ^ c^) G {H{u, Qi) x H{u, Q2)) nL° °(Q ) satisfies 
{4-5)-{4-9) then c defined on Q such that c^^ = c' is i/ie solution of (4.3). In this 
sense the two formulations are equivalent. 



REMARK 1. This analysis can be generalized to the multi-domains case where 

= Uili a, f], n vtj = 0, Q = Uf=i Q^, dQ^ n dQj = r.^, q, ^n,x (o, t), and c e 

iJ(u, Q) is a solution of (4-3) if and only if {c\q^ , cig^: ' ' ' j '^Iq^) — ("^^i ' ' ' i c^) € 



J^j^-j^ _ff (w, Qi) is a multi-domains STILS solution, i.e 



[ Vt • (uc')Vt • {u<f') dxdt = / /Vf • {inp') dxdt i = 1, • • • , P 



= c> 



on 

on f,, ^0, J = !,■■■, P. 



for all ipi G Ho{u, Qi.dQ^ U 



In the sequel is some constant which do not depend on j. 

4.2.3. Iteration-by-subdomain scheme . We propose now a parallel iterative 



procedure to solve the two-domains STILS formulation (4.5|-(4.9|. Consider for i 
1,2 



V- = G Fo(«,Q.,5Q,-), V, • iu^,) G L-(Q,), G i'(r-, |(Sf|n)|da)}, 



G-=7-(^-)- 

For j > 0, assuming (c^'^ , c^'-' ) G G2 x Gj^, find c^'-'+^ and 0^'^+^ respectively 
the least-squares solution of 

(4.16) Vt • (Sc^^J+i) = / 

(4.17) ci'J+i=Cb,i for (a;, i) on aQ]- 

(4.18) ciJ+i=c2.J for(x,t)onf- 
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and 

(4.19) 
(4.20) 
(4.21) 



Vt • (uc- 



2j + l\ 



/ 



c^'^^^ — Cb^2 for (a;, t) on 



for (a;,t) on r+. 

In the STILS sense these subproblems admit unique solutions. To analyze the conver- 
gence of this substructuring algorithm, let us remark that the subproblem ( 4.16 1-( 4.18 1 
can be decomposed into STILS-homogeneous and Dirichlct-homogeneous parts 



^ Ir- 



where the STILS-homogeneous part c^^i+'^'^i g Hq{u, Qi,dQi ) is defined as the least- 
squares solution of 



(4.22) 
(4.23) 
(4.24) 



Vt • (mc^'^'+I'^'^) =0 

^0 for (x,t) ondQ^ 



and the Dirichlet-homogeneous part c^'* € Hq{u, Qi,Ti ) is the least-squares solution 
of 



(4.25) 
(4.26) 
(4.27) 



Cb,l 

c^'* = 



for [x, t) on dQi 
for (x, t) on r~. 



Assuming u £ L'^{0,T; BVinY), Vt • m 0, we have Vt • (uc''/^'^'^ 
(theorem 4.3 1 and 



|Vt • {uc 



1J + 1.A^2^||2 



)lli^(Q.) + ll^ 



< lie 



in L^(Qi 



.2 J" II 2 

r-,(u|n)' 



By the trace inequality (proposition 4.4 1 , the Curved- Poincare inequality (theorem 



4.5) and the proposition 4.7 



T+,(«|n) 



r-,(«|n)' 



with A^(T) = imin(C(T),l) 



In the same way, there exists a unique least-squares solution of (4.19l-(4.21 1 given 
by c^J+i = c2 J+i'M + c5 e H{u, Q2) such that 



r+,(«|n) 



Thus the following estimate result is proved. 

Proposition 4.11. There exists a constant > depending on T such that 



riJ + l||2 



'r+,(u|n) 



<M,i^ 1,2 for all J > 0. 
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Therefore there exists a subsequence of c^'^ (still denoted by this symbol) con- 
verging to some weakly in H{u,Qi). The same argument is applied to the sequence 
c^'-' . This allow us to pass to the limit as j tends to infinity in the STILS forms of 
(4.16 )-(4.21 ). This limit is a solution of the two-domain STILS formulation (4.51- 

Thus the following convergence result is proved. 

Theorem 4.12. Under the assumptions of the proposition the sequence 
(c^'^ j c^'^) defined by the suh structuring algorithm {4^1(w-{4-21) admits a subsequence 
which weakly converges in H{u,Qi) x H{u,Q2) to (c^,c'^) which satisfies the two- 
domain STILS form {4--5)-(4-9)- 



REMARK 2 . If T 



or r+ — 0, the sub structuring iterative scheme (4-.16}- 
(4-21) converges in one iteration. 

5. Time-marching approach. To simplify the notations we suppose cqo = 0, 
the general case being obtained by superposition. Let now decompose the time domain 
(0, T) into M parts 



and put 



in 



< ti < <tM =T 



M 



Q=[jQ\ Q'' ^nx itk-i,tk) 



k=l 



In the space-time description the flow has a constant direction at the inte rface dQ C] 
dQ^ = SI X then with the same assumptions of proposition 4.6 it's seen that 

c e II{u,Q) is the solution of (4.3) if and only if for k = 1,---,M, Ck = ci^^ e 
II{u,Q'') is the unique solution of 



(5.1) 
(5.2) 

for all <fik e HQ{u,Q^,dQ^'-) 



Vt • {uck)Vt ■ (uifk) dtdx 



/V( • (uipk) dtdx 



Ck{x, tk) = Ck-i for X on Q 



In the same way (c\c^) £ II{u,Q^) x II{u,Q^) is the solution of (4.5)-(4.9l if 
and only if for k = 1,---,M, {cl,cl) = (c|q,,cJ^qJ g H{u,Q\) x H{u^) is the 

unique solution of 



(5.3) 

(5.4) 
(5.5) 

(5.6) 
(5.7) 



/ Vt • (u4)Vt • (w(^^) dtdx = / /V( • (w(^^) dtdx 



c\{x^tk) — c\_i for x on ill 



/ • (uc^)Vt • (m^2) ^^^^ ^ f f^^.iu^D^i^^ 
c1{x, tk) — c^-i for X on VL2 



for all e Ho{u,Q\,dQl^-) x i/o(u, 0^, ^g^'") 
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Where 



Q = Oi U Q2, Q, ^n,x {tk-i,tk), i 



1,2, andr'-' = ag^nag^ 



Using theorem 4.10 one shows that the muhi-domains forms (5.1 1-( 5.2 ) and (5.3 )- 
(5.7) are equivalent. Moreover for k — 1,---,M, the sohition of (5.3)-(5.7) can be 
obtained by the iterative scheme (4.16 )-(4.21 ). 

In the sequel we show how to separate the two dimensions space and time to solve 
the above multi-domain problems with finite element discretizations. 

6. Finite element discretizations in time and space. In this section we 
assume that Vf • u = 0, and the equation (3.1 1 becomes 

(Vtc|5) ^ f. 



Then the different formulations (4.3) and (4.5)-(4.9) are now performed by separated 
finite element discretizations in the following way. 

For the space dimension, we consider {Vhi^)}h be a family of regular and quasi- 
uniform triangulations of fi, i.e satisfying the following conditions: there are two pos- 
itive constant a and a independent oi h = max;ce'Ph(o)^K; such that for all /C G Vh(S^) 
(i) h/c > crh 

(a) the angles of JC are bounded from below by a 

where hjc is the diameter of /C. This triangulation can be induced by two inde- 
pendent grids Vih and V2h defined on the subdomains fli and fl2, and which are 
compatible on F (that is they share the same edges therein). 

Let {ipi ■ ■ ■ iPn} be the basis of a finite dimensioned subspace Vh C H{u, Q), where 
tpi is the polynomial function which interpolates the pairs 

(xi, 0), (2:2, 0), • • • , (xi, 1), • • • , {xn, 0), {xi,X2, ■ ■ ■ ,Xi, ■ ■ ■ , xpf are the nodes of the grid 

For the time dimension, the domain [0, T] is divided into n subintervals 
[t°,t^], [t^,t^], as follows 

and for each {t^~^,t''), k — 1, 



T 



T 

with t & [t''~'^,t% T ^ t'' - t''-\ 

C[t^~^,t''] generated by a^^i and a^. 



• , n we consider affine basis functions 

ak-iit) = -{t''-t), ak{t) = -{t-t''~^) 



Let C2 be the 2— dimensional subspace of 



Put now 



= rj, X i = l,2 



Then the TMQl (Time-Mar chin g Method with Ql basis functions in time) approx- 
imation [6] of the problem (4.3) consists in finding for > 1, Ch = c'^ & 
V ® C2) n H{u, Q^") such that 

(6.1) / tCh\u)(W t^Phlu) dtdx = 
Jqi>' 

(6.2) Chix,t''-') = c'^-' 



f{Vt^h\u) dtdx 



for X onD, 
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for all e i?''-" = (V* (g> Ci) n i7o(w, Q, dQ^" -). 

Consider {V'l''' , ■ • • , V'^i } the set of the basis functions associated with the nodes 
lying on ili, i = 1,2. And be the finite dimensioned subspace induced by this set 
of functions. Then the TMQl approximation of (4.5)-(4.9) is expressed as follows: 
for k > 1, (4,4) = {c1,4) e fff X H'i = ((Vf ® C2) n h{u,qC)) x ((V^^ ® C2) n 



Q2 )) such that 



(6.3) 

(6.4) 
(6.5) 

(6.6) 
(6.7) 



Ik 



(Vt4|2)(Vt^i|u)dtdx= / /(Vt^i|S) dtdx 



c\^^ for X on Vti 



(Vt4|2)(Vtv?fJu) dtdx = / ^ fi^M^u) dtdx 



for , 



^2 



{ipl ^l) e i/f'" X = ((vf ® C2) n i/o(^i, aQ( u r^'"-)) x ((v^' ® c^) n 

where '~ and '+ are respectively the inflow and outflow part of F^ . 
For each time subdomain (t'^^^,^'^), the equivalent result between the discrete prob- 
lems (6.1)-(6.2) and (6.3)-(6.7) is given by 

Theorem 6.1 . As suming V t ■ u — and f — in (2.1), the n if C h e is 
a solution of (6.1)-(6.2) then {ch\Q^ , ^hi pj € (-ff? x H2 ) satisfies (6.3)-(6.1). Re- 
ciprocally if (cjj, 4) € X Hk satisfies (6.3)-(6.1) then Ch defined on Q'' such that 
Ch\Qik = is a solution of { 6.1)-{ 6.2). 

Finally as in the continuous level (section 5) the above TMQl approximations 
are equivalent. 



In order to solve the two-domains discrete formulation (6.3)-(6.7), it's convenient 
to apply the iteration-by-subdomain scheme. 



Thus for j > until convergence, find cj^ 
such that 



1 j+i 



c1^'+' e and cl'^^ 



(6.8) / {Vtc];''+'\u){VtA\u) = / ^ fi-^tvWu), for all ^1 G 77^'° 



(6.9) 
(6.10) 

and 

(6.11) 

(6.12) 
(6.13) 



Ir/fc,- 



^'^'{x,t^-') 



2,i+i 



2j- 



fiyt^l\u), for all ^ i7. 
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Following the same guideline of the continuous level, we get the estimates 



<X(t), i= 1,2, 



and the discrete convergence result is given by 

Theorem 6.2. Under the assumptions of the p r oposi tion the sequence 
(cjj'-' , c^'"* ) defined by the sub structuring algorithm (6.8)-[6.13) admits a subsequence 
which weakly converges in x H2 to {c\, c^J which satisfies the discrete two-domains 
STILSform (K^-jtl^. 



Higher orders time basis functions can be considered in the above analysis. For 
example one can to apply the time-marching approach with (52— basis functions [B] 
(TMQ2) which are defined on [0, r] by 



afe-i(i) = 



afe(i) 



(t- 


m- 


-r) 


(0- 


i)(o 


-r) 


{t- 


0){t- 


-r) 


(i- 


m 


-r) 


(t- 


0){t- 


-i) 


(r- 







REMARK 3. The subproblems (6.8)-(6.10) and (6.11)-(6.13] can be approximated 



differently both in space and time. Indeed: 

1. For the space domain, the two grids Vih andV2h can be choosed independently. 

2. For the time subdomain {t^~^^t^), the choice of the basis functions can also be 
different for the two subproblems, for example one can to adopt the TMQl method for 
{6.8)-{ 'KT^ , and the TMQ2 approximation for (6.11)-(6.13). 



In the next section some numerical results are presented to analyze the potential 



of the iterative scheme (|6.8|)-( |6.13| ) to solve the STILS-DDM formulation (|6.3[)-(|6.7[). 
In this sense the STILS-DDM numerical solution and the STILS numerical solution 
are compared for the cylinder rotating example |21] . 

7. Numerical results. We consider the TMQl approach to analyze numerically 
the STILS and STILS-DDM formulations 



7.1. Time integration. To solve (6.1 1-(6.2) by a TMQl method, we choose tph 
independent of afe_i. So iph cancels on t — t^^^ . This is motivated by con side ring 
the local problem in the domain 

n X {t''-^}. If 



n X with an initial condition (6.2) on 



c''{x,t) 



N 

J2Mx){ct~^ak-i{t) + c'[ak{t)), 
1=1 



problem (6.1) becomes in this case, after time integration see [B] 

1 



(7.1) 



N 

1=1 

N 



1 



1 



dx 



1=1 



-(Vi/;,|u)(VV',|u) 



1 



Z Z T 



dx 
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for alH in {1, • • • , iV}. For the STILS-DDM algorithm ([6!8|)-( |6.13| , we choose ip]^, i 
1,2 independent of a^-i, and if 



TV" 



1=1 



then (6.8)-(6.10) and (6.11 1-(6.13) become 



N" 



(7-2) Y.^S" 



1=1 



-(vv^) '\u)tpi +--0;' yi 

I T 



1=1 



m = 1, 2, for all i in {1, • • • , N}, and j > 0. 



. 6 
dx. 



-(v^rV)(vvrV) 



The related sparse, symmetric, and positive define systems (7.11, (7.2) are solved 



by means of the Preconditioned Conjugate Gradient (PCG) algorithm. 

7.2. Cylinder with a slot example. Let consider the cylinder with a slot 
example taken from [21\ The domain is the square =] — 1, Ip and discretized in 
100 X 100 elements, of type Ql. The initial condition is 



c(a:,y,0) 



1 if (|x| > 0.05 or y > 0.7) and i? < 3 
elsewhere 



where R = \/x^ + {y — 0.5)^, the velocity field has the form 

v{x,y,t) = {-y,x), 



and the final time is T = 27r. 

At final time, i.e. after a complete rotation, the exact solution coincides with the 
initial condition. The domain is decomposed by two subdomains 

H - 1, o[ X ] - 1, i[, =]o, i[ X ] - 1, i[, r = {o}x] - 1, i[. 



Here are the results for numerical solutions of equation ( |2.1| using the STILS and 
STILS-DDM methods. Thus for r = 27r/800, h = 2/100, the figures [tI] fL2\ and 
7.3 show respectively the snapshot of the numerical solutions at 100th, AOOth, 500th 
step on 800. The cylinder rotating effectuates perfectly. In more the shape of the two 
slots are similar, and they are preserved compared to the initial condition. Thus this 
example of simulation justifies well that the time-marching approch of the STILS- 
DDM procedure even its simplest form (Ql time basis functions ) is well encouraged 
to approximate the STILS problem. 
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STILS - DDM 





Fig. 7.2, 



TMQl numerical solutions of STILS (left), STILS-DDM (right) at AOOth step 
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STTLS STILS - DDM 




Fig. 7.3. TMQl numerical solutions of STILS (left), STILS-DDM (right) at 500th step 



STILS STILS - DDM 




Fig. 7.4. TMQl numerical solutions of STILS (left), STILS-DDM (right) at final time 

REMARK 4. The STILS-DDM solution presents some small oscillations on the 
interface which doesn't appear in the STILS numerical solution see fig. \7.S\ These 
numerical oscillations are induced by the solving of the two subproblems to obtain the 
interface condition. However with a small tolerance these oscillations are not so suf- 
ficient to create a major difference between the two schemes. 

The two methods present also some differences and difficulties induced by the 
two finite element discretizations. The following results show the influence of the 
discretization parameters on the different algorithms. 

7.2.1. Influence of the discretization parameter r. The relative variation 
which is the maximal variation of the integral concentration across the time steps, 
divided by the minimal value measures in an inverse way the global conservativity. 
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Ideally it should be since there is no input nor output of concentration and the exact 
solution inside the unit disk is simply a rotation of the initial condition. 



T 


Number of iter 


(min,max) 




(min,max) 


relat variat 


relat variat 




< 


STILS-DDM 




STILS 


STILS-DDM 


STILS 


27r/50 


2 


(-0.33 , 0.91) 


(- 


-0.35 , 0.92) 


3.14- 10-^ 


2.65-10-2 


27r/100 


3 


(-0.39 , 1.06) 


(- 


-0.41 , 1.07) 


1.60- 10-2 


3.98-10-3 


27r/200 


3 


(-0.23 , 1.25) 




-0.28 , 1.25) 


4.98- 10-3 


8.70-10-* 


27r/400 


4 


(-0.28 , 1.28) 




-0.19 , 1.28) 


1.97-10-3 


3.59- lO-'' 


27r/800 


4 


(-0.23 , 1.21) 




-0.19 , 1.20) 


7.44 • lO-'' 


1.44 - 10-4 



Table 7.1 
Numerical results respecting on r 



0.035^ 




0.00 0.05 0.1 



Fig. 7.5. Conservativity respecting on t 



STILS - DDM. T = 2k/400 



STILS - DDM, T = 2x/2(X) 
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STILS - DDM. T = 2k/ 100 



STILS - DDM, T= 2k/S0 





Fig. 7.6. STILS-DDM-TMQl numerical solutions at final time for 



values of t. 



As showed on the table 7.1 and fig. 7.5 considering a fixed value oih — 2/100 the 
global conservativity for each scheme decreases (relative variation increases) when r 
increases, with a light advantage of STILS. In the STILS-DDM case this is illustrated 



by a deterioration of the shape of the slot at final time (see fig. 7.6) 



7.2.2. Influence of the discretization parameter h. The same remarks are 
observed when h takes enough high values and r = 27r/800, then the global conserva- 
tivity decreases (see table 7.2 and 7.3, fig. 7.7 ) 



Number of 


Number of 


Numb of iterations 


(min,max) 


relative variation 


elements 


nodes 


< 


STILS-DDM 


STILS-DDM 


1600 


1681 


3 


(-0.25,1.20) 


1.2 10"^ 


2500 


2601 


3 


(-0.20,1.30) 


0.89 10-3 


3600 


3721 


4 


(-0.18,1.26) 


1.17 IQ-^ 


10000 


10201 


4 


(-0.23,1.21) 


7.44 10-'' 


40000 


40401 


5 


(-0.19,1.23) 


4.06 10-5 



Table 7.2 

Numerical results of STILS-DDM-TMQl respecting on h 



Number of 
elements 



Number of 
nodes 



(min,max) 
STILS 



relative variation 
STILS 



2.76 10-^ 
5.99 10-3 
0.98 10"^ 
1.44 10"^ 
1.87 10-5 



1600 1681 (-0.33,1.21) 

2500 2601 (-0.17,1.28) 

3600 3721 (-0.21,1.23) 

10000 10201 (-0.19,1.20) 

40000 I 40401 I (-0.15,1.20) 

Table 7.3 

Numerical results of STILS-TMQl respecting on h 
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0.030^ 




0.010 0.015 0.020 0.025 0.030 0.035 0.040 0.045 0.050 



Fig. 7.7. Conservativity respecting on h 

Finally one remarks also that the two schemes induce some small overshootings 
and undershootings which are nearly the same in the two cases. This due to the fact 
that the TM scheme is not conservative nor even TVD. At end it can be to remark 
that the maximal number of iterations for the STILS-DDM algorithm to obtain con- 
vergence varies very little with t or h (see table 7.1, and 7.2). 

These numerical results lead to the following remark 

REMARK 5. Experimentally, in oder to product a good quality of the solution, it's 
necessary to have < ^. It's look like an unstable condition while the time-marching 
approach of STILS is unconditionally stable see fl8\/ . 
When > \ there is a phenomenon look like a numerical diffusion. 

8. Conclusion. Using the STILS method a new space-time domain decompo- 
sition algorithm to solve scalar conservation laws is presented. A convenient way to 
analyze numerically this algorithm is to discretized separately by finite element meth- 
ods the space and time dimensions. The different grid parameters are very influent 
on the convergence of the iterative scheme. The Hansbo example illustrates well the 
results of this work i.e to approximate the STILS problem by a domain decomposition 
method. Higher orders and all sort of combinations with the techniques in Galerkin 
variants are possible for the proposed iteration-by-subdomain scheme. 
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